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We report recent progress of non-perturbative formulation of supersymmetric Yang-Mills. Al- 
though lattice formulations of two-dimensional theories which are fine tuning free to all order in 
perturbation theory are known for almost ten years, however, there were only few evidence for 
the validity at non-perturbative level. In this talk we argue that most numerical studies so far do 
not capture the physics in continuum, and add new evidence that lattice formulation works at non- 
perturbatively. We further point out that, by combining two-dimensional lattice and matrix model 
techniques inspired by D-brane dynamics in superstring theory, a non-perturbative formulation of 
the four-dimensional maximally supersymmetric Yang-Mills theory, which is free from the fine 
tuning at least to all order in perturbation theory, is obtained. 
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1. Introduction 

Supersymmetric Yang-Mills (SYM) theories play prominent roles in theoretical particle physics. 
In particular, maximally supersymmetric theories are expected to describe the non-perturbative na- 
ture of superstring/M theory [[T], ^ ^, |5]]. Lattice formulations of SYM theories are pursued be- 
cause many important question would be answered only via numerical simulations. However, it is 
not a straightforward task, because the supersymmetry cannot be preserved completely on lattice: 
the SUSY algebra contains infinitesimal translation, which is broken on lattice by construction. 
Therefore, even if a given lattice theory converges to a supersymmetric theory at tree level, SUSY 
breaking corrections may arise radiatively. In order to control the divergence one needs some exact 
symmetries at discretized level. In 4d ^ = 1 pure SYM, by keeping the chiral symmetry one can 
obtain the correct supersymmetric continuum limit [^'. In several extended SYM theories, it is 
possible to keep a few supersymmetries unbroken. In two dimensions, those exact supersymme- 
tries, together with other global symmetries, forbid SUSY breaking radiative corrections at least 
to all order in perturbation theory [^, ^ [T^, [TT]]. (A similar statement holds for the Wess-Zumino 
model, and numerically tested in [p^.) 

In one dimension (i.e. supersymmetric matrix quantum mechanics), the situation is much eas- 
ier. Because the theory is UV finite, one does not have to rely on exact symmetries and hence a 



simple momentum cutoff prescription works [13]. In fact, as demonstrated in [13], the momentum 



cutoff method is much more powerful than a usual lattice regularization and detailed Monte-Carlo 
studies have been done. In particular, the gauge/gravity duality between DO-brane quantum me- 
chanics and type IIA superstring theory [||] has been tested, and beautiful agreement including the 
a' correction has been confirmed [ |T4| , [TSf]^. Furthermore long-disrance physics relevant for the 
matrix theory conjecture has also been studied [|T6|]. 

In order to study higher dimensions, one must establish fine-tuning free formulation of those 
theories. Obvious first step is non-perturbative test of two-dimensional theories. In particular one 
should see whether the fine-tuning-free nature persists to non-perturbative level. In § |2| we show 
numerical results which strongly support the validity of these formulations at non-perturbative 
level. Based on the success of two-dimensional lattice formulation, in § ^ we provide a hybrid 
formulation of two-dimensional lattice and matrix model techniques, which enables us to put four- 
dimensional maximally supersymmetric Yang-Mills theory on computer. 

2. Non-perturbative test of two-dimensional lattice SYM 

Numerical simulations of the supersymmetric gauge theories often suffer from the fermion 
sign problem^ . One of the exceptions is 4d ^ = 1 SYM and its dimensional reduction. Here we 
study 2d .yV = (2,2) SYM, which is the dimensional reduction to two dimensions. The action is 

S = j[^dx j^dy Tr^^Fl, + \{D^Xif -\[XuXjf - ly/r^'D^i/.- ^v/r[X,-, | , (2.1) 



'For recent numerical studies, see [^. 

^Other numerical studies in the context of the string theory can be found in ||l7[[r^, [l9| ], for example. 
^In the case of maximally supersymmetric matrix quantum mechanics, agreement with the dual gravity prescription 
has been observed by ignoring the phase of the Pfaffian, even when the sign fluctuates violently [jl^. 
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where jU and v run x and y, i and j run 1 and 2, and = (F^ , F' ) are gamma matrices in four 
dimensions. X,- are N x N hermitian matrices, are x A/^ fermionic matrices with a Majorana 
index a and the co variant derivative is given by = — i[A^, ■ ]. The only parameters of 
the model are the size of circles L,- and Ly. (Note that the coupling constant can be absorbed by 
redefining the fields and coordinates. Therefore we take the 't Hooft coupling A to be 1. Then the 
strong coupling corresponds to the large volume.) 

One obstacle for the simulation is the existence of the fiat direction, along which two scalar 
fields Xi and X2 commute. In contrary to a theory on M^'^, there is no superselection of the moduli 
parameter in this case. That is, eigenvalues of scalars are determined dynamically. Therefore, some 
mechanism which restrict eigenvalues to a finite distribution is necessary for the stable simulation. 
In addition, to obtain an interesting dynamical system, having a (small) finite region for the eigen- 
values is important as well; if the eigenvalues of the scalar spread so large, the theory would run 
into the abelian phase, which is just a free theory"^. In this work, we introduce soft SUSY-breaking 
mass to scalar fields 

H'^N j (f-x £ TrXf, (2.2) 

!=1,2 

so that the fiat direction is lifted. It is crucial to control the fiat direction for various reasons. We 
have just mentioned two of them — stability of the simulation and interesting non-abelian phase. 

In the Weyl notation, with an appropriate choice of the gamma matrices, the Dirac operator is 
written as 

D = io^D^, (2.3) 

where a" = — /I2 and a'(/ = 1,2,3) are Pauli matrices. By using a^(/a'^)a^ = (ia^)* and the fact 
that is real in adjoint representation, we obtain 

o^Do^=D*. (2.4) 

Therefore, if (p is an eigenvector corresponding to an eigenvalue A, o^(p* is also an eigenvector, 
with eigenvalue A*. They are linearly independent and eigenvalues appear in a pair (A, A*). This 
assures the positivity of the determinant after removing A = modes. 

At discretized level, positivity of the determinant can be lost. In Cohen-Kaplan-Katz-Unsal 
(CKKU) model [^] an existence of the sign problem had been reported. However it turned out 
that these studies did not capture the continuum physics, because the lattices were too coarse. In 



[20] and [21], CKKU model and Sugino model [10] have been studied and it has been confirmed 



that the determinant becomes real positive in the continuum limit (see also [22]). In Fig. |I| we 
have plotted the distribution of the phase of the pfaffian^ in the SU (3) Sugino model. Here we 
chose Lx = Ly = L and imposed periodic boundary condition for all fields along both directions. 
In the left panel we fixed the physical volume to be L = 0.707 and used three lattice sizes, 4x4, 
5x5 and 6 x 6. It turns out the peak around zero becomes sharper as the lattice becomes finer. In 
the right panel we have fixed the lattice size and changed the lattice spacing. It can be seen that 



^Which phase is preferred is in fact a dynamical question. At large-A', tlie flat direction is lifted and the system stays 
non-abelian phase This phase is an analogue of the black 1-brane solution in type IIB supergravity. 

^Here we adopt the Majorana fermion. Hence we calculate the pfaffian rather than the determinant. 
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Figure 1: [Sugino] Argument of the Pfaffian in SU{T>) theory. The scalar mass is jU = 0.20. The left panel 
is for a fixed volume 0.707 x 0.707 and thus different lattice spacings. The right panel is for a fixed 4x4 
lattice with various physical volumes (thus various lattice spacings). (Figure from [|T]]) 




Figure 2: The expectation value of the Wilson loop at /i = 1.0. Extrapolation to the continuum limit 
has been performed. The gauge grope is t/(2) and SU (2), respectively. (Figure from \V^\^ 



smaller lattice spacing gives sharper peak. Both plots clearly show that the sign disappears as the 



continuum limit is taken. For more details, see [21 1 



Since we have established the absence of the sign problem in 2d = (2,2) lattice SYM, we 
can perform various tests. So far, conservation of the supercurrent in the Sugino model has been 
observed in [23]. Here we add another result. In Fig. § we have plotted the expectation values 



of the Wilson loop winding on compactified circle, which are calculated in CKKU and Sugino 
models. Two models give the identical result, which strongly suggests that they converge to the 
same continuum limit without performing parameter fine tuning. 

Although theories with more supersymmetries suffer from the sign problem, it is a technical 
problem attached to the importance sampling and we can expect the restoration of all supersymme- 
tries in the continuum limit also for them. Further investigation along this Une is desirable. 



3. Four-dimensional SYM out of two-dimensional lattice 

In four dimensions, a few exact supersymmetries are not strong enough to control the radiative 
correction, although it does reduce the number of fine tuning parameters down to three [^]. In 



[ 25 ] it has been pointed out that the fine tuning problem can be circumvented by using the matrix 
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model techniques^ . The basic idea is simple. Consider the Berenstein-Maldacena-Nastase matrix 
model 1^], which is a deformation of the DO-brane matrix quantum mechanics keeping maximal 
supersymmetry. It has a fuzzy sphere solution, which is BPS. If one considers /:-coincident fuzzy 
sphere, three-dimensional U {k) non-commutative super Yang-Mills (NCSYM) on fuzzy sphere is 
obtained. By taking appropriate \wge.-N limit, NCSYM on flat non-commutative plane is obtained, 
and by turning off the non-commutativity, one arrives at usual SYM on [p9|]. This method can 
easily be generalized to four-dimensional theory; if one starts with 'BMN-like' two-dimensional 
theory with fuzzy sphere solutions, one can construct 4d theory. The crucial point is that such 
two-dimensional theory can be regularized by using the lattice. Bosonic part of the action of this 
two-dimensional theory is given by [ p^ ] 

S2d = -^ f d^x Tr\ ^2 + ^ ^D^Xif - j [Xj,Xjf + ^ £ + 2i^X, [X2,X,] 

-^X6[X7,X8]}. (3.1) 



The deformation by /i breaks 14 out of 16 SUSY softly. As shown in p5[], this model can be 
regularized by lattice keeping two supercharges unbroken, and the continuum theory is obtained 
without parameter fine tuning, at least to all order in perturbation theory. 



The continuum action (3.1) has constant BPS fuzzy sphere solution^ 



Xa{x) = ^La (fl = 1,2,3), Xi{x)=0 (/ = 4,.--,8), (3.2) 
where L„ are M x M matrices satisfying SU (2) commutation relation 

[l-a-,l^h] = i^ahcLc- (3.3) 

By taking ^-coincident fuzzy sphere solution. La = L^^^^^ Ijt, where l}a^^^^ is the {M /k) x (M/k) 
irreducible representation, we obtain 4d U (k) theory on fuzzy sphere. Essentially, adjoint action of 
La is identified with the derivative and [Xa, ■ ] is regarded as the gauge covariant derivative [^]. 
The noncommutativity is given by ~ k/{iJL^M) and UV/IR momentum cutoffs along spherical 
directions are jiM/k and /i, respectively. 4d coupling is given by = 27i6g\^. In order to get 
continuum 4d theory, we take large-M and small pL limit while fixing k and g\^. In that limit, 
maximal supersymmetry is restored because soft SUSY breaking parameter pL goes to zero. One 
can take a limit with any value of non-commutativity Q, and — > limit is expected to be smooth 
in the maximally supersymmetric theory. That the limit should be smooth is natural physically, 
because a possible obstacle is a new IR divergence arising due to the UV/IR mixing reflecting the 
UV divergence, which should be absent in UV finite theories. 

In the above we assumed the radius of the fuzzy sphere does not deviate from classical value. 
Whether it is the case or not should be tested by numerics. If the radius is renormalized, we should 
take it into account by replacing the parameters in the mapping rule with renormalized ones. 



''For another matrix model approach which works at large-W, see 
^This background preserves exact supersymmetries at discretized level. 
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4. Conclusion 



We reported recent progress of the non-perturbative formulations of SYM. Lattice formu- 
lations of two-dimensional ^ = (2,2) SYM has been tested non-perturbatively and confirmed 
to work without parameter fine tuning beyond the perturbative level. It is also shown that four- 
dimensional SYM can be constructed by combining two-dimensional lattice and a matrix model 
method. 

There are various problems one can attack based on above results. Firstly, 2d = (2, 2) 
SYM can have interesting physics on its own (see e.g. [22, 30, 31, 32, Rl]) and further numerical 



study is desirable. 2d maximal SYM is even more interesting because it describes [ 34 ] the black 
hole/black string phase transition p^. It is important to establish the absence of fine tuning in 



this case and perform a large-scale simulation along the line of [17|. Numerical study of four- 
dimensional SYM would be the most important. The formulation explained in § |3| would be useful 
to understand the AdS/CFT duality further, especially away from the weak and strong coupling. 



Another interesting class of theories are SUSY QCD. At large-A'^, the large-A'^ reduction [ ]27[ ] can be 
applied to SUSY QCD [^. It would provide us with a tool to study important problems e.g. the 
study of spontaneous SUSY breakdown or the test of Seiberg duality. It is interesting especially 
because it may have a close connection to physics at LHC. 
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